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Abstract 


Wave  propagation  and  random  media  are  defined  and  the  nature  of  the 
mathematical  problems  arising  in  wave  propagation  in  random  media  is  described. 
The  two  principal  types  of  methods  for  solving  these  problems  -  honest  and 
dishonest  methods  -  are  explained.  These  methods  are  first  illustrated  by- 
considering  the  geometrical  optics  of  a  random  medium  by  one  method  of  each 
type.   Some  new  results  are  obtained  by  an  honest  method  and  some  errors  in 
a  previo\is  work  are  pointed  out.   Comparison  is  made  between  the  results  of 
the  two  methods  and  the  reasons  why  they  disagree  are  explained.  As  a  sec- 
ond illustration  of  an  honest  method,  an  analysis  of  the  reduced  wave  equa- 
tion in  a  random  medi\;mi  is  presented.  Some  known  results  are  obtained  in  a 
new  way  which  is  simpler  than  the  usual  one  and  which  appears  to  be  capable 
of  yielding  further  results. 
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1.   Introduction 

Wave  propagation  is  one  of  the  means  by  vhlch  energy  travels.   'Propa- 
gation' is  the  process  whereby  the  energy  moves  from  one  region  of  space  to 
an  adjacent  region^  and  'wave'  is  a  general  term  for  a  moving  spatial  dis- 
tribution of  energy.   The  matter  in  the  region  of  space  through  which  the 
propagation  occxirs  is  called  the  transmission  medium.   A  random  or  stochas- 
tic medi\jm  is  a  family  of  media  together  with  a  probability  distribution 
which  specifies  the  probabilities  of  the  various  members  of  the  family. 
Thus  wave  propagation  in  a  random  medium  refers  to  propagation  in  each  mem- 
ber of  the  family  of  medla^  together  with  the  probability  of  each  member. 
This  probability,,  when  associated  with  the  wave  motion  in  each  medium, 
characterizes  a  random  wave  motion. 

Mathematically  a  wave  motion  is  described  by  a  vector-valued  function 
u(x,  t)  of  the  position  vector  x  and  the  time  t.  As  a  consequence  of  the 
physical  laws  governing  the  wave  motion,  the  function  u(x,t)  satisfies  cer- 
tain equations.  Usually  they  are  partial  differential  equations  of  hyper- 
bolic type  and  often  of  symmetric  hyperbolic  type.   The  transmission  medixom 
is  characterized  by  a  vector-valued  function  n(x,t)  which  enter  the  coeffi- 
cients of  the  equations.  A  wave  propagation  problem  Is  that  of  determining 
a  solution  u(x, t)  of  the  equations  which  satisfies  certain  auxiliary  condi- 
tions.  These  conditions  are  usually  initial  and  boundary  conditions.   The 
problem  is  said  to  be  well  set,  well  posed  or  properly  posed  if  it  possesses 
a  unique  solution  which  depends  continuously,  in  an  appropriate  nonn,  upon 
the  data  given  in  the  a\ixlllary  conditions  and  upon  the  coefficients  in  the 
equation. 

A  random  medium  is  characterized  by  a  family  of  functions  n(x, t,(jj)  de- 
pending upon  a  parameter  ui.   Here  w  Is  a  point  in  a  probability  space  Q   in 


which  a  probability  measure  dP(u))  is  defined.   Of  coiirse  the  function 
n(x;,t,w)  must  be  measurable  with  respect  to  dP(u)).   If  a  well-posed  wave 
propagation  problem  is  formulated  for  the  random  medium^  it  has  a  solu- 
tion u(x^t,co)  for  each  w.   This  family  of  solutions  describes  a  random  wave 
motion.   The  equations  satisfied  by  u  with  the  coefficient  n(x,t^u))  are 
called  stochastic  equations  -  stochastic  differential  equations  if  they  are 
differential  equations.   Thus  in  this  case  wave  propagation  in  a  randum  me- 
dium is  part  of  the  subject  of  stochastic  differential  equations. 

The  reasons  for  studying  random  wave  propagation  are  exactly  the  same 
as  those  for  studying  any  other  random  phenomenon.   There  are  essentially 
two  such  reasons.   In  the  first  place^  we  may  wish  to  consider  a  case  in 
which  the  medium  is  not  known  precisely^  but  in  which  it  is  known  to  be  a 
member  of  a  certain  family  of  media.   If  we  also  know  the  probability  that 
it  is  any  particular  member  of  the  family^  we  can  determine  the  probability 
that  the  wave  motion  is  any  one  of  the  associated  family  of  wave  motions. 
Then  we  caii  determine  the  expected  or  mean  wave  motion^  as  well  as  its 
variance  and  other  statistics.   This  statistical  information  can  be  used 
to  estimate  what  is  likely  to  be  the  wave  motion  in  the  case  under  consider- 
ation.  In  particular  some  statistical  properties  of  u  may  depend  only  upon 
known  statistics  of  n.   When  we  consider  majiy  media^  all  of  which  are  mem- 
bers of  the  same  family,  we  may  expect  the  statistics  of  the  corresponding 
wave  motions  to  correspond  to  those  of  the  random  wave  motion  occurring  in 
the  random  medium. 

In  the  second  place  we  may  wish  to  consider  a  known  medium  with  very 
complex  properties.   Then  the  associated  wave  motion  will  also  be  very  com- 
plex and  its  precise  determination  may  be  impractical.   Certain  statistical 
properties  of  the  wave  motion  in  a  random  medium,  of  which  the  medium  in 
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question  is  a  member^  may  be  less  complex  and  more  easily  determinable.   In 
other  vords,  It  may  be  easier  to  solve  the  random  problem.   If  the  random 
medium  is  appropriately  chosen,  these  statistical  properties  may  be  closely 
related  to  the  actual  properties  of  the  vave  motion  in  the  known  medium 
under  consideration. 

We  shall  see  that  there  are  two  kinds  of  methods  used  to  solve  problems 
of  wave  propagation  in  random  media.   These  may  be  called  'honest'  and  'dis- 
honest' methods  respectively.   In  an  honest  method  the  solution  u(x,t,co)  is 
first  determined  for  each  value  of  w.   This  solution  may  sometimes  be  found 
exactly  and  explicitly,  but  more  often  it  is  expressed  in  the  form  of  a 
series  in  some  parameter,  or  as  a  sequence  of  iterates,  or  by  some  other 
approximation  procedure.   In  the  process  of  solving  for  u(x, t,co)  randomness 
plays  no  role  and  therefore  it  proAAides  no  advantage.   The  second  step  in 
an  honest  method  is  to  compute  the  mean  value  of  u(x,t,u)),  as  well  as  its 
variance  and  other  statistics,  from  the  explicit  expression.   In  this  step 
randomness  may  have  the  helpful  effect  of  yielding  simpler  expressions  for 
the  statistics  of  u  than  those  for  u  itself. 

In  a  dishonest  method  randomness  is  utilized  before  u(x,t,uj)  is  deter- 
mined. For  example,  if  the  mean  value  <u>  is  sought,  the  original  equations 

for  u  may  be  'averaged'  to  yield  equations  for  <u>.   In  such  cases  the  aver- 

2 
aged  equations  also  involve  terms  like  <u  >  or  <nu>.   Dishonesty  enters 

2       2 
throxigh  the  assvmrption  that  <u  >  =  <u>  or  <nu>  =  <n>Cu>  or  some  similar 

unproved  assumption.   The  reason  for  making  such  assijmptlons  is  to  obtain 
determinate  equations  for  <u>.   Other  examples  of  dishonest  methods  will  be 
described  later.   In  all  cases  probability  is  introduced  before  u  is  deter- 
mined and  an  \mproved  assumption  is  made  about  some  statistical  property  of 
the  random  wave  motion.   The  assimiption  usually  simplifies  the  problem  so 


that  It  becomes  solvable. 

Dishonest  methods  have  the  advantage  over  honest  ones  that  they  simplify 
the  problem  to  be  solved.  As  a  consequence,  a  problem  vhich  can  be  solved 
honestly  only  by  a  perturbation  method  might  be  solved  dishonestly  without 
the  use  of  a  perturbation  expansion.   Then  the  dishonest  solution  would  be 
applicable  for  all  values  of  the  relevant  parameter  while  the  honest  solution 
would  be  valid  only  for  small  values  of  it.   For  this  and  similar  reasons, 
many  of  the  significant  and  non-trivial  results  in  the  theory  of  wave  propa- 
gation in  random  media  have  been  obtained  by  dishonest  methods.  Many  of 
these  results  have  compared  well  with  experiment.   Thus  one  of  the  important 
mathematical  problems  in  this  field  is  to  justify  the  dishonest  methods  by 
showing  that  their  results,  in  some  sense,  are  approximations  to  honest  solu- 
tions. A  clear  vmderstanding  of  the  circumstances  in  which  this  is  the  case 
would  permit  the  introduction  of  many  more  useful  dishonest  procedures. 

Most  of  the  work  on  this  subject  has  been  done  since  19^5^  having  been 
stimulated  by  practical  problems  of  radio  wave  propagation  through  the  atmos- 
phere and  ionosphere,  sound  wave  propagation  in  the  ocean  and  the  atmosphere, 
light  transmission  through  the  atmosphere,  etc.   The  recent  book  by  L.A. 
Chernov'--'  contains  a  rathe'r  complete  bibliography  and  an  understandable 

account  of  the  present  state  of  this  subject.  Additional  material  is  con- 

[21 
tained  in  the  book  by  V.I.  Tatarski  •- -^  .   The  related  subject  of  random  wave 

motion  in  a  non-random  medium  is  surveyed  by  M.  Bom  and  E.  Wolf  '-  ■*        -^ . 

Some  other  work  is  contained  in  the  symposium  volumes  edited  by  W.C.  Hoffman'--' 

and  Z.  KopalJ-  J  .   Because  of  the  complete  bibliographies  in  these  books,  we 

shall  give  relatively  few  references. 


Physicists  have  studied  extensively  the  propagation  of  waves  through 
random  collections  of  discrete  scatterers.   A  clear  formulation  of  a  prob- 


lem of  this  type,  together  vlth  a  new  dishonest  method  for  treating  it,  was 
given  by  L.L.  Foldy^  -1  .   This  method  wr.s  generalized  to  other  problems  by 
M.  Lax'-  -' ,  who  has  also  reviewed  much  of  the  previous  work  in  this  field. 

Many  problems  have  been  treated  by  this  method  and  others  by  V.  Twersky'-  -' . 

[91 
Recently,  J.  Bazer'-  -'  proved  that  for  a  one-dimensional  scattering  problem, 

this  method  does  yield  the  correct  result.   I.  Kay  and  R.A.  Silverman'-  -' 
investigated  the  accuracy  of  another  method,  the  single  scattering  approxi- 
mation, by  determining  the  extent  to  which  randomness  reduces  the  importance 
of  multiple  scattering. 

The  next  section  describes  how  wave  propagation  in  a  random  medium 
applies  to  the  problem  of  the  twixikling  of  a  star.   In  Section  J  we  consl-. 
der  light  rays  in  a  slightly  inhomogeneous  medium,  and  in  Section  ^  we  apply 
these  results  to  a  random  medium  obtaining  some  new  results.  We  then  compare 
these  honest  results  with  corresponding  dishonest  ones  in  Section  5,  where 
we  treat  a  light  ray  in  a  random  medium  as  a  Mai^off  process.   In  Sections 
6  and  7  we  use  the  results  of  Sections  3  sixd  k   to  determine  the  phase  and 
amplitude  fluctuations  of  a  wave  in  a  random  medium.   Our  results  are  exten- 
sions of  previously  known  ones.   This  completes  our  discussion  of  geometri- 
cal optics  in  a  random  medium  .   In  the  final  section  we  present  a  brief  new 
treatment  of  the  reduced  wave  equation  in  a  random  medium.   This  section  is 
independent  of  the  preceding  ones. 

2.   An  application 

Before  Illustrating  the  techniques  used  to  analyze  wave  propagation  in 
random  media,  we  shall  describe  a  physical  phenomenon  in  which  such  propaga- 
tion plays  a  role.   It  is  the  scintillation  or  twinkling  of  a  star.   The  most 
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appropriate  theory  of  vave  propagation  to  describe  this  phenomenon  is  geo- 
metrical optics.   According  to  this  theory^  light  travels  along  certain 
straight  or  curved  paths  called  rays.   These  rays  are  determined  by  ordinary 
differential  equations  in  which  appears  the  index  of  refraction  n(x, t),  a 
scalar  function  vhlch  characterizes  the  transmission  medium.   This  medium  is 
the  earth's  atmosphere  in  the  present  case.  If  n=l  the  rays  are  straight  lines 
which  emanate  from  the  star  in  all  directions.   One  of  them  (really  a  narrow 
beam)  enters  the  eye  of  an  observer  who  is  viewing  the  star.   The  direction 
from  which  the  ray  enters  the  eye  is  the  apparent  direction  of  the  star. 
Actually  n(x^t)  differs  from  unity  by  a  small  amount.  As  a  consequence,  the 
rays  deviate  slightly  from  straight  lines  and  enter  the  observer's  eye  from 
slightly  different  directions  at  successive  Instants  of  time.   Therefore  the 
star  appears  to  be  moving  about  its  mean  position.   Its  apparent  Intensity 
also  fluctuates. 

We  would  like  to  calculate  the  apparent  direction  and  intensity  as  func- 
tions of  time.   To  do  so  it  would  be  necessary  to  know  the  index  n(x,t).   This 
is  practically  impossible  becatise  the  variation  of  n(x_,t)  with  time  and  posi- 
tion results  from  the  trubulent  motion  of  the  atmosphere.   Consequently  analy- 
sis would  appear  to  be  impossible.   In  the  face  of  this  difficulty  we  treat 
the  atmosphere  as  a  random  medium.   The  random  medivmi  must  be  so  chosen  that 
its  Important  statistical  properties  correspond  to  measurable  properties  of 
the  atmosphere.   In  tMs  way  we  give  up  the  possibility  of  calculating  the 
apparent  direction  and  intensity  at  any  particular  time.   Instead  we  can  calcu- 
late statistical  properties  of  the  apparent  direction  and  intensity  which  may 
be  related  to  the  actual  temporal  distribution  of  apparent  directions  and  in- 
tensities.  We  shall  consider  this  example  further  in  the  next  section.   In 
doing  so  we  shall  make  use  of  the  fact  that  a  ray  traverses  the  atmosphere  so 
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qiilckly  that  the  index  does  not  change  significantly  during  the  traversal. 
Consequently  we  may  assvune  that  the  index  is  independent  of  time  in  deter- 
mining the  rays. 

5.   Light  rays  in  a  slightly  inhomogeneous  medium 

Let  n(x, e)  denote  the  index  of  refraction  of  a  random  medium,  vhlch  may 
be  written  as 

n(x,e)  =  1  +  e  |a(x)  .  (l) 

The  quantity  €  measures  the  deviation  of  the  index  from  unity.   The  index  also 
depends  upon  a  parameter  co  which  we  shall  not  write  explicitly.  We  wish  to 
determine  the  ray  x(s,€)  which  starts  from  the  origin  in  the  direction  of  the 
unit  vector  u.   Here  s  denotes  arclength  along  the  ray.   The  ray  also  depends 
upon  w.   The  equations  which  x  satisfies  are 

{n^')'    =  Vn  (2) 

x(0)  =  0  (5) 

x'(0)  =  u  ,  (u^  =  1)  .  (4) 

These  equations  have  a  unique  solution  which  depends  continuously  upon  n  and  u 
in  an  appropriate  norm,  if  n  is  continuously  differentiable.  Thus  the  problem 
(2)  -  (k)    is  well  posed. 

To  find  X  we  shall  determine  its  derivatives  with  respect  to  €  at  e  =  0 
and  then  express  x  by  means  of  its  Taylor  series  in  e.   Thus  we  shall  employ 
an  honest  method,  in  the  first  phase  of  which  probability  plays  no  role.  Let 
us  first  set  €  =  0  in  (2)  -  (^4-)  and  denote  x(s,0)  by  x  (s).   Then  we  obtain 

x"  =  0  (5) 

o 

x^(o)  =  0  (6) 

x;(o)  =  u  .  (7) 


The  solution  of  (5)  -  (7)  is 
X  (s)  =  us  . 


(8) 


Nov  ve  differentiate  (2)  -  (k)   vith  respect  to  e,    set  e  =  0,  and  denote 


X  (s^O)  by  X  (s)_,  obtaining 


x"  =  V  ^(x  ) 

€         '  O' 


x'  •  V  ^(x  ) 

0     ^  O 


(9) 


X  (0)  =  x'(0)  =  0 


(10) 


The  right  side  of  (9)  is  just  the  component  of  Vp.  vhich  is  normal  to  x'  =  u. 
Let  us  call  it  the  transverse  gradient  and  denote  it  by  V  p .  Then  the  solu- 
tion of  (9)  and  (10)  is 

s 

x^(s)  =  J  (s-t)  V^  n(ut)dt.  (11) 

o 

Differentiating  (2)  -  {k)   tvlce  with  respect  to  e  at  e  =  0  yields 

x-(s)  =  2(x^  .  V^)V^  ^^(x^)-V^  H^(x^)-2^;x^  .  V^  H(x^)-2^^x;  •  V  n(xj   (12) 


X  (0)  =  x'  (0)  =  0 


The  solution  of  (l2)  and  (15)  is 


(15) 


x^^(s)  =  2  /  (s-t) 


""e  '  ^T  ^  ■  I  "^T  ^^-^^^T  ^^^-^^^'   '^t^)^T  ^ 


dt.     (l4) 


In  the  integrand,  the  argument  of  |i  is  ut.   Thus  to  the  second  order  in  e  "we 

have 

^        2  ^ 
x(s,e)  =  x^(s)  +  e  x^(s)  +  ^  \^i^)   +  0(€^)  .  (15) 
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k.       Light  rays  In  a  random  medium 

Let  us  now  take  account  of  the  randomness  and  compute  some  statistical 
properties  of  x(s,€).  Let  us  begin  with  the  mean  value  <x(s,€)>  which  is 
just  the  s\am  of  the  mean  values  of  the  terms  on  the  right  side  of  (15)  • 
Since  x  (s)  =  us  is  independent  of  n  and  therefore  of  w,  <x  (s)>  =  us.   From 
(11)^  by  interchanging  the  order  of  taking  the  mean  value  with  integration 
and  differentiation,  which  we  assume  to  be  permissible,  we  obtain 


<x 

e' 


Xs)>=  /"(s-t)  V^  <^(ut)>dt.  (16) 

o 

In  (l6)  only  <n>  occurs.  We  see  that  it  is  not  yet  necessary  to  specify  the 

probability  space  n,  the  measure  dP((jo)  nor  n(x,w)  if  instead  we  give  <ii(x)>. 

We  therefore  assume  that 

<(i(x)>  =  0  .  (17) 

From  (l6)  and  (l?)  it  follows  that 

<5c  (s)>  =  0  .  (i8) 

To  compute  the  mean  value  of  x  (s),  which  is  given  by  (l^),  we  must 
know  the  mean  values  of  certain  quadratic  expressions  in  ^  and  its  deriva- 
tives.  They  can  be  determined  from  the  correlation  function  <(j(x  )|i(x  )> 
which  will  also  be  reqxiired  in  later  calculations.   Therefore  we  assiune  that 
it  is  given  by 

<n(x^)|i(x2)>  =  <^^>  N(  Ix^-x^l )  .  (19) 

2 
Here  the  mean  square  fluctuation  <n  >  is  assumed  to  be  a  constant  and  the 

correlation  coefficient  N  is  assumed  to  be  a  function  of  the  distance  |x^-Xp| 

only.   By  the  Schwartz  inequality  it  follows  that  |n|  SI.   Of  course  N(0)=1. 
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We  aisc  require  that  N'(0)  =  0  and  that  N(c»)  =  0.   It  follovs  that  N  fO)  ^  n. 

rr  ' 

In  addition  N  should  become  small  beyond  a  distance  a  vhich  ve  shall  call  the 

correlation  distance. 

2  2       ? 

From  the  assumption  that  <u  >  Is  constant  ve  have  <V  u  >  =  V  <m  >  =  0. 

Thus  the  mean  of  the  second  term  In  the  Integrand  In  (ik)   is  zero.   To  com- 
pute the  mean  of  the  first  term  ve  use  (ll)  for  x  and  the  first  term  becomes 

s  '^2  r  ^ 

I J     (s-tgjCtg-t^)  J[y^^n(ut^)  .V^2  V^^I^Cutg)  V  dt^dt^  .  (20) 


o 


The  expression  in  brackets  can  be  revritten  as  (V   •  V  )V  M(utp)n(ut  )  and 

since  V  ,  =  -  V  ^,  vhen  axipllei  to  functions  of  t^-  t.   its  mean  value  is 
Tl      T2-^  2   1 

If  ve  set  hg-t^l  =  r  then  N  =  N(r)  and  V^^  =  ^m-   Since  V^  N(r)  is  a  func- 
tion of  r  onlv^  its  transverse  gradient  is  zero.   Therefore  (21)  and  the  mean 
of  (20)  are  zero.   In  a  similar  vay  ve  find  that  the  mean  of  the  last  term  in 
(l^)  is  also  zero.   Thus  only  the  third  term  remains  and  ve  obtain 

s 

<^(s,e)>  =  '^  -ue^  Ms-tXv^  ^)S  dt  +  0(€^)  (22) 

o 

2  2 

To  evaluate  <(V     |i)    >  ve  have,    since  V       N  =  -  V  p  N  and  since  V    N(0)   = 

|v^(o), 

<(V^  m)  ^>  =  <n^>       11m       V      .  V       N(  I X  -x^  I ) 

=  -  </>  V^^N(O)    =  -  I  </>  V^N(O)    .  (23) 

2 
Thus   finally^    taking  account  of  the   fact  that  V    N(0)   =  5N     (O)^   v  e   obtain 

^(s,6)>  =  uQ  +  €^  <mS  N^^(O)s^)  +  0(€^)    .  {2k) 
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From  (24)  we  can  solve  for  the  length  s  such  that  a  ray  of  this  length 
travels  a  mean  distance  L  in  its  original  direction.   This  value  of  s  is 

s  =  L  -  e^  <^^>  N  (O)L^  +  0(€^)  . 
rr^        ^ 


(25) 


Equation  (2^)  gives  the  mean  location  of  the  endpoint  of  a  ray  of  length 

-♦ 
s  which  starts  from  the  origin  in  the  direction  u.   This  mean  position  is  on 

the  straight  line  through  the  origin  in  this  direction,  which  we  expect  to  be 
the  case  from  symmetry.  However,  its  distance  from  the  origin  is  less  than  s, 
as  we  also  expect.   The  terms  given  in  (24)  do  not  suffice  to  compute  <x(s,€)> 
for  very  large  values  of  es.   This  can  be  seen  by  noting  that  for  such  values, 
the  mean  position  of  the  endpoint  would  lie  in  the  direction  -u,  according  to 
(24)  since  N  (O)  S  0.   This  resiilt  is  unreasonable,  as  is  the  fact  that  the 
mean  distance  from  the  origin  to  the  endpoint  would  exceed  s  for  very  large 
€S.  We  conclude  that  the  approximation  (24)  for  <5(s,€)>  is  not  uniform  in  s. 
Th\is  more  terms  in  the  series  are  reqiilred  for  large  values  of  es.   They  in- 
volve third  and  higher  order  correlation  functions  and  moments  of  |j. 

Let  us  now  compute  the  mean  value  of  x(s,€)    .   From  (I5)  we  have 


-i2 


x(s,e) 


2    2 

=  S   +  € 


X  (s) 


+  e  su  •  X  (s)  +  0(e)  . 
ec       ^   ' 


(26) 


In  writing  (26)  we  made  use  of  the  fact  that  u  •  x  =0.   By  \iBing  (ll)  and 
(l4)  in  (26)  we  obtain 


-i2 


x(s,e) 


2    2/ 
=  s  +  e 


(s-t)V^ndt 


-Se^s 


J  {s-t)iV^p.fd.t   +  O(e^) 


(27) 


The  mean  value  of  the  third  term  on  the  right  side  in  (27)  has  been  computed 
in  reducing  (25)  to  (24).   To  compute  the  mean  of  the  second  term  we  write  it 
as  follows 
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s  s 


<^>  J  J    (s-t^)(s-t2)V^^  .  V^2  N(  |ut^-  utgl  )dt^dt. 


o  o 


8  S 


^^^  //^^■*1^^^"^2)  "^T^  Ndt^-t^Ddt^dtg  . 


o  o 


(28) 


We  now  define  r  and  r  by 

o  -^ 


t^  =  r^  +  r/2     tg  =  r^  -  r/2  . 
Then  (28)  becomes^  if  we  define  N(-r)  =  N(r)^ 

o  s+(r/2) 


<(x  f>=   -  </> 


-s    -r 


/2 


/    ^2  r 


V^  N(r)dr  dr 
T  ^    o 


(29) 


2 
<U  > 


s  s-(r/2)r 
o    r/2 


"■o' 


2  r 


V^  N(r)dr^dr  .    (30) 


Upon  performing  the  integration  with  respect  to  r  and  combining  the  resulting 
Integrals  we  obtain  from  (50) 


<(x^)2>  = 


<U^>  f    v|  N(r) 


2   2s' 
rs  +  ^j- 
5 


dr 


(51) 


To  simplify  (31)  we  note  that 


v|N(r)  =  2r"^  W^(r) 


(52) 
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With  the  aid  of  (52),   ve   can  rewrite   (51)   as 


2      /-    x2 
e     <(x^)    > 


2        2 
e     <|i  > 


N^(r) 


2r 


o  2       i^s- 
2s     +  ^— 
3r 


dr 


(55) 


Now  (27)   becomes 


1—    s 


,-»,        ,i2  2  2        2 

<|x(s,e) I    >  =  s     -   e     <^  > 


/    \(r) 


^f--  Ss^H-i^^dr-  23'b^^(0) 


»-0 


+    0(€^) 


(5i^) 


For  values   of  s  large   compared  to  the   correlation  distance  a,    (5^)    can  be  sim- 
plified to 


,--         ,,2  2225 

<|x(s,e)  I    >  Si  s  -   e     <n  >  s 


-     /     r"-"-  N  (r)dr  -   2N     (O) 
5   J  r^  rr^ 


+  0(€O..(55) 


Equation  (3^)  or  (55)  gives  the  mean  square  distance  between  the  end- 
points  of  a  ray  of  length  s.  Let  P  denote  the  transverse  displacement  of  the 
endpolnt  of  a  ray  of  length  s  from  the  straight  line  tangent  to  the  ray  at  its 

initial  point.   Then  equation  (55)  gives,  to  order  e  ,  the  mean  square  value 

2 
<p  >.   For  s  large  compared  to  the  correlation  distance  a,  (55)  yields  for 


<.=> 


2 

<p  >  ;s. 


1  2   2   5 
4-e  <n  >  s-^ 


go 


N  dr  +  0(€"') 


(56) 


By  Bubstraction  of  (56)  from  (55)  we  obtain  the  mean  square  distance  which 
the  ray  travels  along  its  original  direction.   For  s  large  compared  to  a,  it 
is 

<|u.x(s,€)|^>  Sis^  +  2e^  <U^>N^^(0)s^  +  O(e^)  .  (57) 


-> 


Ik  - 


Let  us  finally  calculate  the  mean  square  value  of  x'(s^e)  -  u.  This  is 
the  mean  square  value  of  2  -  2  cos  a(s^e)  where  a  is  the  angle  between  x'  and 
u.   From  (15)  we  have 

(x'  -Zf   =  e^{x^f  +  0{e^)    .  (58) 

From  (58)    and   (ll)   we  obtain 


<(x'    -  uf>  =  e^  <  (      /  V^  n  dt  )      >  +  O(e^) 


(39) 


By  proceeding  as  In  (28)  -  (33)  we  find 


<(x'  -  u)^>=  2M  -  <cosa(s,e)>j  =  he^   <n^> 


For  s  large  compared  to  a  this  becomes 


N(s)-l-s  /  r'hi   dr 


+  0(e 
(^0) 


3n 


<(x'  -  u)  >=  2(  1  -  <cosa(s,€)> 


he^   <H^>  s  /  r'-Sr^dr  +  O(e^)  . 
o 

(^1) 


If  the  right  side  of  (^l)  is  small  compared  to  imlty^  we  have 

2 
2(1  -  <cos  a>)  ^  <a  >.   Then  (^l)  gives  the  mean  square  value  of  a. 

The  results  (36)  and  (4l)  can  be  written  simply  in  terms  of  the  ray  dif- 
fusion coefficient  D  which  is  defined  by 


D  =  -  e^  <[i^>     /  r"-^  N^(r)dr 


Then  {%)    and  (^l)  become 

2    i|.    5      ^ 
<p   >  ^a^   Ds-^  +  0(g). 

2  3 

<a  >  ^  4-Ds  +  0(6  ). 


(42) 


(^5) 
(kk) 
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An  incorrect  form  of  (^1)  vas  derived  by  L.A.  Chernov  [l^p.lY^eq. Jl] • 
His  method  of  derivation  is  the  same  as  that  given  here  but  he  did  not  intro- 
duce a  small  parameter  e  and  expand  systematically  in  povers  of  it.   As  a 
consequence  he  made  an  error  in  going  from  his  eq.  (25)  to  eq.  (26)  on  page 
16.   Therefore  his  result  eq.  (51)  and  his  diffusion  coefficient  eq.  (5^) 
page  17  are  incorrect. 

5.   Light  rays  as  a  Markoff  process 

In  the  preceding  section  we  treated  light  rays  in  a  slightly  inhomogene- 
ous  random  medium  by  an  honest  method.  We  computed  the  mean  values  of  some 
quantities  associated  with  the  rays  and  obtained  simple  expressions  for  them 
when  the  ray  length  s  was  large  compared  to  the  correlation  distance  a.   Al- 
though the  results  are  correct,  they  are  not  useful  for  very  large  values  of 
s,  as  can  be  seen  by  examining  them.   We  must  conclude  that  the  quantities  we 
have  calculated  in  powers  of  €  are  not  uniform  in  s. 

To  obtain  results  which  are  valid  for  large  values  of  s  a  dishonest  method 
has  been  employed L J .   it  consists  in  treating  the  ray  direction  x' (s)  as  a 
random  function  of  the  arclength  s.   The  fundamental  unproved  assumption  on 
which  the  method  is  based  Is  that  x'(s)  is  a  continuous  Markoff  process. 
This  means  that  the  probability  that  the  ray  ha^  a  given  direction  at  s  +  ds 
depends  only  upon  the  probability  distribution  at  s.   Let  us  denote  this  pro- 
bability density  by  P(9,^s).   Then  as  a  consequence  of  conservation  of  pro- 
bability, P  satisfies  the  Chapman-Kolmogoroff  equation.   The  differential 
form  of  this  equation  is  the  Fokker-Planck  equation 


P  = 
s 


^   (sin  9  P,),  +   ^ 


sin  9  '       9'9    ,  2^ 

sin  9 


'  li-   ^  •  (^5) 

As  ->-o 


2  _^ 

Here  <h,  >  is  the  mean  square  value  of  the  angle  a  between  x'(s)  and  x'(s+As) 
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when  x'(s)  is  given.   We  have  calculated  it  in  the  preceding  section  and  it 

2  2 

is  given  by  (^O) .  From  (^O)  it  follows  at  once  that  <a  >  is  of  order  (As  ) 

for  small  As  and  thus  the  limit  in  (k-^)   is  zero.   However^  if  we  use  (kk)   for 

2  ^ 

<a  >    instead  of  (4o),  and  neglect  0(e  ),  the  limit  in  (^5)  is  D  and  (^5)  be- 


comes 


sin  9 


(^6) 


It  would  seem  to  be  inconsistent  to  use  (^^),  which  is  valid  for  s 
large  compared  to  the  correlation  distance  a,    and  then  to  let  s  tend  to  zero. 
But  this  kind  of  procedure  is  common  in  applied  mathematics.   It  can  usually 
be  Justified  by  interpreting  it  as  determining  an  asymptotic  expansion  with 
respect  to  some  appropriate  parameter.   In  the  present  case  this  is  so  and 
the  parameter  is  the  correlation  distance  a.   To  show  this  we  shall  compute 
the  limit  of  s~  a<ix  >  as  a  tends  to  zero^  by  means  of  (^O).   In  doing  so  we 
first  observe  that,  from  dimensional  considerations,  N  is  a  function  of  a  r, 
say  N(r)  =  M(a""'"r).  Now  (4o)  yields 


-1    2  2   2   -1 

lim  s  a  <a,  >  =  11m  ^e  <|i  >  s   a 


a-»o 


ar-MD 


M(a  s)-l-s 


/   r  a  m' (a  r)dr 


2   2-1 
=  ^e  <n  >  s   lim 


arK) 


-1 
a  s 


aM(a  a)-a-s 


j      x"^'(x) 


dx 


ke^  <^%    I     x"-S{'(x)dx 


i^l) 


—  1    2 
From  (^7)  and  {h2)   we  see  that  for  any  s  7^  0,  (^s)"  <x  >  -  D  as  a  tends  to 

zero.   Then  the  limit  in  (^5)  is  D.  We  now  see  that  (^6)  is  valid  in  the 
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limit  in  which  the  correlation  distance  a  tends  to  zero.   Thus  the  solution 
of  (^6)  is  the  leading  term  in  the  expansion  of  P  for  small  values  of  a. 

We  have  seen  that  the  assumption  that  x'(s)  is  a  continuous  Markoff  pro- 
cess leads  to  (^5).  When  a  tends  to  zero  (^5)  becomes  (^6).   If  P  is  inde- 
pendent of  ^,  which  is  the  case  for  all  s  if  it  is  so  for  s  =  0,  (^6)  further 
simplifies  to 

''s  ■  iis^  ^^'- «  V,  •  ■  («) 

Let  us  now  solve  (^8)  subject  to  the  'initial'  condition 

P(0,9,^)  =  5(9)  .  (kg) 

This  'initial'  condition  corresponds  to  the  ray  at  s  =  0  certainly  pointing 
in  the  direction  9  =  0.   The  solution  of  (^8)  and  {k-<^)   is  expressible  in  terms 
of  Legendre  polynomials  P  (cos  0)  as 


Hs,Qj)   =1^  2  (2n+l)Pjcos  9)e-^^^-^l^°^  .  (50) 


From  (50)  we  see  that  P(<»,9,(^)  =  l/k-n   so  that  for  large  s  all  directions  be- 
come equally  probable.   From  (50)  we  easily  find 

-2Ds  ,   , 

<cos  9>  =  e     .  (51) 

Let  us  compare  (51)  with  the  result  for  <cos  a>  which  can  be  obtained 
from  (^0).   We  observe  that  a  and  9  are  defined  in  the  same  way  so  that  the 
two  results  should  agree.   However,  they  do  not  agree.  Since  (ko)   was  ob- 
tained by  an  honest  method  and  (5I)  by  a  dishonest  one,  we  might  be  tempted  to 
conclude  that  the  dishonest  method  has  yielded  an  incorrect  result.   But  the 
preceding  discussion  in  this  section  shows  that  this  is  not  the  proper  conclu- 
sion. Instead  we  must  recognize  that  the  two  results  have  different  domains 
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of  validity.  The  result  (k-Q)  is  valid  for  small  values  of  s  while  (51)  is 
valid  for  large  values  of  s.  If  we  expand  (kl)  for  s  large  and  (5I)  for  s 
small  we  obtain  agreement  to  two  terms : 

<cos  a>  :^  1  -  2Ds  ~  <cos  9>  .  (52) 

This  agreement  shows  that  the  two  domains  of  validity  overlap  slightly.  Of 
course  the  question  of  whether  the  result  (5I)  is  valid  at  all,  since  it  is 
based  on  the  Markoffian  assumption;,  is  not  answered  by  these  considerations. 

The  result  (5I)  can  be  used  to  determine  the  mean  square  distance  be- 
tween the  endpoints  of  a  ray  of  length  s.   To  use  it  for  this  purpose  we 

write 

s  s 

<|x(s)l^>  =  <  J   x'(t^)dt^-  J    x'(t2)dt2  > 


s  s 


<x'(t^)-x'(t2)>  dt^dtg 


o  o 


s  s 


=   /  /  <cos  9(t^,t2)>  dt^dtg 


o  o 


s  s 


=  J  J    expC-SDlt^-t^D  dt^dt^ 
o  o 

=  D  -  -2  C'  -  -    J    '  (55) 


For  small  s  this  becomes 
n2 


< 


x(s) 


P   2   ^ 
>  5i  s  -  -  Ds-^  .  (5^) 
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The  first  term  of  this  result  does  agree,  but  the  second  does  not  agree, 

x(s)   >  for 


vlth  (55)  >  which  gives  the  previously  calculated  vsilue  of  < 


large  s.   That  result  involves  N  (O)  in  addition  to  D,  while  all  the  re- 
sults of  the  present  method  involve  only  D.   We  must  conclude  that  the 
domains  of  validity  of  (55)  and.  (5^)  do  not  overlap. 

Let  us  finally  compute  </x(s)-uJ  >  by  the  present  method.   This  is 
the  mean  square  value  of  the  distance  vhich  the  ray  travels  along  its  ini- 

tiail  direction.   It  is  given  by 

s  s 
<(x(s).u)^>=   /  /  <cos  e(t^)cos  9(t2)>  dt^dtg  .  (55) 

o  o 
By  using  (50)  we  finally  obtain 

:i;  s^  -  2Ds^  .  ■'  (56) 

The  second  form,  of  (56)  holds  for  small  Ds.   Its  first  term  agrees  with 
(57);  which  gives  the  same  quantity  computed  by  the  previous  method.   How- 
ever, the  second  term  does  not  agree. 

2 
By  substraction  of  (56)  from  (55)  we  obtain  <p  >,  the  mean  square  trans- 
verse displacement  of  the  endpoint  of  the  ray  from  its  initial  tangent  line. 
It  is  given  by 

Si  ^  Ds^  .  (57) 

The  second  form  of  (57),  which  holds  for  small  Ds,  agrees  with  (^5)^  obtained 
by  the  previoios  method. 
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These  calculations  suffice  to  illustrate  this  method  and  to  shov  to 
vhat  extent  its  results  agree,  or  disagree,  with  those  of  the  preceding 
honest  method. 

6.   Phase  fluctuations 

In  geometrical  optics  rays  do  not  occur  individually  "but  only  in  normal 
congruences.   A  normal  congruence  is  a  tvo  parameter  family  of  rays,  all  of 
vhlch  are  normal  to  some  smooth  surface.   Every  surface  normal  to  the  rays  of 
such  a  congruence  is  called  a  vavefront.  With  every  normal  congruence  of 
rays  we  associate  a  complex  valued  field  v(x)  which  we  write  as 

v^)  =  A(x)  e^^*^^^).  (58) 

The  possibly  complex  factor  A(x)  Is  called  the  amplitude  of  the  field  and 
(^(x)  is  called  its  phase.   These  quantities  are  uniquely  defined  in  terms  of 
V  by  the  requirements  that  <!>   be  independent  of  the  propagation  constant  k 
and  that  A  depend  upon  k  at  most  through  a  factor  of  a  power  of  k.   The  pro- 
pagation constant  k  is  equal  to  2it/X  where  k   is  the  wavelength  of  the  field 
in  a  region  of  refractive  index  equal  to  unity.   The  amplitude  A(x)  may  be 
either  a  scalar  or  a  vector.   The  real  time  harmonic  field  ■V(x,t)  determined 


,-*■.    -ivt 
v(x)e 


Here  t  denotes  the  time  and  v 


by  v(x)  is  given  by  V(x,t)  =  Re 
the  angular  frequency  of  the  field. 

Let  X  =  z(a,  (3)  be  the  equation  of  a  wavefront  S^,  i.e.  a  surface  ortho- 
gonal to  a  normal  congruence  of  rays.  We  shall  suppose  that  the  parameters 
a  and  p  are  so  chosen  that  at  a  , p  the  following  equations  hold 


o  o 


,2   ,   -  -   „  -    -1  -*   -    -1 


(^a)^=  ^^p)^=  ^'  vv  °^  V  ^r  V  V  '2^  'p'  '-^'p  ^  "•      ^^^^ 
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Here  u  denotes  the  unit  normal  to  S  while  p  and  p  are  the  principle  radii 

of  S  at  a  , 3  .   These  equations  imply  that  at  a  ,  3  the  parameters  a  and  3 

measure  arclength  along  the  lines  of  principle  curvature  on  S  . 

We  denote  by  x(a, 3,s)  the  ray  normal  to  S  at  a, 3'   Thus  measuring  s 

from  S  .  ve  have 
o' 


x(a,3^0)  =  z(a,3) 
x'(a,3^0)  =  u(a,3) 


(60) 
(61) 


The  ray  through  a  given  point  y  is  determined  by  the  values  of  a,  3  and  s  such 
that 

y  =  x(a,3.s)  .  (62) 

Once  these  values  are  determined,  we  define  the  phase  ^(y)  by 


My)  =  J   n 


x(a,3^t) 


s 

dt  =  S  +  €    /   [X 

o 


x(a,3,t) 


dt 


(65) 


We  have  taken  6   to  be  zero  on  S  • 

^  o 

To  evaluate  ^  we  must  first  solve  (62)  which  we  shall  do  by  expansion  in 
powers  of  €•   By  using  (15)  we  can  write  (62)  as 


y  =  x^(a,3^s)  +  €  x^(a,3,s)  +  0(€  ). 


Here  x  (a, 3,s)  is  given  by 


X  =  z(a, 3)  +  su(a, 3) 


(6^) 


(65) 


Let  us  set  e  =  0  in  (6^+)  and  denote  the  solution  by  a  ,3  ,s  .   Then  by  dif- 
ferentiating (6^)  with  respect  to  e  and  setting  e  =  0  we  obtain 

0=ax  (a.3,s  )+3x„(a,3^s  )+sx'(a,3.E  )+x(a,3,s  )  . 
e  oa^  o'^o'  o"^   *^e  03^  o'  o'  o'    c  o^  o"^o^  o     e^  o'^o'  o 


(66) 
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By  using  (65)    and  (59)    snd  the   fact  that  u.x     =   0  we  readily  obtain  from 
(66)   the  results 


a     = 
e 


-z   'X   (a  ,  3   ;S    ) 
a     e^    00     o 


1  +  s 


-1 
0^1 


(68) 


1  +  s^p-1 


(69) 


s     =  0 


(69) 


We  now  use  these  results  in  (63).   Since  i-i  =  V^.(x  +  a  x  +  3  x  „)  for 

^      '  6         6   e  oa   €  03 

e  =  0^  we  obtain 


X   (a  ,3  ,t) 
o  00 


dt 


o 


x  (a  ,3  ,t) 
o^  00 


z  -x  (s  ) 
X  (t)  -  ^L-^° 

1  -^  %pl 


-  (1  +  tp-^)z^ 


TI  (l  +  tp^  )z 

1  -^  ^^2  - 


dt  +  0(e'^) 


(70) 


We  next  -use  (ll)  for  x  in  (70)  and  compute  the  mean  value  of  ^(y).  A 
lengthy  but  rather  straightforward  calculation  yields 


2    -=0 


<^(y)>=  s^+  e^  <^i^>(p  ^ 


(71) 


If  p     and  p     are  both  finite  then^    for  large  s^^    (7I)   becomes 


<i!^(y)>  a«  s   (1  +  €  <M  >)   +  0(e'^) 


(72) 
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Thus  the  mean  phase  increases  linearly  with  s     as  if  the  medivim  had  the  con- 

2     2 
stant  index  of  refraction  n*  =  1  +  e  <n  >. 


The  variance  ct 


^G) 


is  given  by 


^(y) 


=  < 


-i2 


2 


> 


2    -.2 
o   ^e   2  '^ee    o   '^e   2  '^ee 


=  <  |s„+  €^g+  ^  j!^^^-  <s^+  e^^+  -^  (^^^>  >  +  0(€^) 


2  ,2       "5 
=  €  <p  >  +  O(e^) 


(75) 


From  (70)  and  (75)  we  have 


^G) 


s   s 
o  o 


-'  J  J  ^^ 

o  o 


-o^^) 


xjt^) 


>  dt^dtg  +  0(6'') 


=  2e^  <4i^>     /  (s^-r)N(r)dr  +  O(e^) 


(7*^) 


For  large  s  this  becomes 
^       o 


^G) 


2  2 
^  2€  <ii  >  s 


/  N(r)dr 


(75) 


In  concluding  our  discussion  of  the  phase  ^,   ve  should  observe  that  ^  is 
a  solution  of  the  eiconal  equation 

(V^)^  =  n^  .  (76) 


It  is  a  solution  which  satisfies  the  boundary  condition 


^  =  0     on  S 


(77) 


There  are  two  such  solutions  which  merely  differ  in  sign.   Thus  the  analysis 
of  this  section  has  concerned  the  solution  of  a  boundary  value  problem  for  a 
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random  partial  differential  equation  of  first  order.  The  analysis  was  based 
upon  the  solution  of  the  equation  by  means  of  the  characteristics,  which  are 
the  rays  in  the  present  case. 

7.   Amplitude  fluctuations 

The  amplitude  A(x)  also  satisfies  a  first  order  partial  differential 
equation  in  which  j^(x)  occurs.   For  a  scalar  A  it  is 

V.(A^(!^)  =  0  .  (78) 

Let  us  seek  the  solution  which  has  the  prescribed  value  A  (a, p)  on  the  sur- 
face  S  given  by  x  =  z(a, 3).   Thus 


z(a,p) 


=  A°(a,p)  .  (79) 


These  equations  can  be  solved  for  A  by  the  method  of  characteristics.   Because 
of 'the  special  form  of  (78)  the  solution  can  be  obtained  easily  by  integrating 
(78)  over  a  volume.   Now  the  volume  integral  of  the  divergence  of  a  vector  is 
equal  to  the  surface  integral  of  the  normal  component  of  the  vector.   To  sim- 
plify this  surface  integral  we  choose  as  the  volimie  a  region  bounded  by  a  tube 
of  rays,  the  wa^/efront  S  and  another  wavefront  S.   Since  V^  is  tangent  to  the 
rays,  the  surface  integral  over  the  tube  vanishes.   Thus  only  the  integrals 
over  the  two  portions  of  the  wavefronts  remain  and  we  obtain 


/ 


A^|Vi!$|dcr-   /  A^|ViZ^|dcr^  =  0  .  (80) 


S'  S' 

o 

In  (80)  S'  and  S'  are  the  two  portions  of  the  wavefronts,  da  and  dc  are  the 

elements  of  area,  and  we  have  made  use  of  the  fact  that  V^  is  normal  to  the 

wavefronts,  pointing  Inward  on  S  and  outward  on  S.   Since  (80)  holds  for  every 

choice  of  S',  and  the  corresponding  S'  determined  by  the  rays,  we  can  conclude 

that 


25   - 


A-'lV^I    =  a''{0)\V^(0)\ 


do- 
o 

dcr 


(81) 


Here  dcr  /dc  is  the   ratio  of  corresponding  area  elements   on  S      and  S,    and  the 
argtiment   zero  denotes   s  =  0  which  corresponds   to  a  point  on  S    .      Since 
\V^\    =  n,   ve  obtain  from  (8l)    and   (79)   the   result 

y n|x(a,3,s)]  y  \      / 


(82) 


To  evaluate  A  we  must  compute  dtf/dcr  ,   which  is  the  Jacobian  of  the  map- 
ping of  S  onto  S  by  means  of  rays.   This  aapping  is  given  by  the  function 
x(a,  p,s)  with  s  having  the  constant  value  corresponding  to  S<.   If  a  and  3 
are  defined  as  before  at  a  ^ p  on  S  ^  then  along  the  ray  through  this  point 
we  have 


do- 
dCT 


=  X  (a  ,3  ,s)  X  x„(a  ^3  ,s)   . 
(_a^  0^0       3^  0*^0   J 


(85) 


If  we  use   (6^)    for  x(a, 3^s),    (83)   becomes 


dor 
do- 


[(l-fp-^s)z^+  ex  J   X    [(l+p-\)l  +  e\~] 


+   0(6^) 


(1+P^  s)(l+P2  s)u  +  €[(H-P^  s)(z^xx^    )   +  (1+P2^s)(x^^x  z   )] 


=  (l+p^^s)(l+P2^s)   +  €[(1+P^^s)u.(z^x  x^    )   +  (l+Pg  s)u.(x^^xz   )]+0(e    ) 


=   (l+P^^s)(l+p"^s)   +  €[(l+P^^s)Zp.x^p+  (l+P2^s)z^.xJ-fO(e^)    .      (Ok) 
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Now  (82)  yields  for  A  the  result 


x(a,P,s) 


=A^(a,3,s)  ^1- 


(l+p-^s)-^lp.x^p+(l+p-^)-^1^.x 


€a 


-  n(z)+n(x) 


+  0(€^). 


(85) 


Here  we  have  introduced  A  which  is  defined  by 

o 


A^(a,p,s)  =  A°(a,|3)(l+P^^s)-^/2(l+p-^B)-^/2 


(86) 


To  obtain  A  at  a  fixed  point,  say  y,  we  must  insert  for  a, 3  and  s  in 

2  2 

(85)  the  solutions  of  (62).   They  are  a  +  €a  +  0(a  ),    p  +  ep  +  0(e  )  and 

2 
s  +  0(e  ).   Then  (85)  yields 


A(y)  =  Aja^,3^,sJ  {1   -  ^ 


(I+P2  s)   Zp.x^p+  (1+p^  s)  z^-x 


ea 


-  m(z)+h(x) 


+  eA 


A  a  +  A  „p 
oa  e   op  € 


+  0(6'").  (87) 


From  (87)  we  find  at  once  that 


<A(y)>=  A  (a  ,3  ,s  )  +  0(e' 
^■' '     o  o  o  o     ^ 


(88) 


Next  we  compute  the  variance  of  A(y) . 


2       2 
<A  >  -  <A>  =  < 


A(y)  -  A^ 


>+0(€^) 


(89) 


We  shall  evaluate  it  in  the  special  case  In  which  A  =  A  ^=  0.   Then  after  a 

^  oa   op 


lengthy  calculation,  we  find  for  s  »  a. 
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2     2  ^   5.2 


12 


1      ^_i_)     ,2(-::i-+     ^2 


5     V  P^+s        pg+s 


P-L+s        Pg+s 


+  s 


P^+Pg 


2  ^  2 

_(p-L+s)  (p^+s)(p2+s)         (Pg+s) 


Here 


(90) 


C  =  J   j"    (v|)^(r)dr  =    J     \r 


r     N     -   r"\ 
rr  r 


dr   . 


(91) 


When 


Pj_  =  Pg  =  P^    (90)   simplifies  to 


2     2        3  2 

o  o     €  <M  x;s  A 

<A^>  -  <A>^  =  


3(p+s) 


0        7    s  ^2 

2—       ■(    —  +  sp  +   2p    ' 


(92) 


If  the  initial  wavefront  is  a  plane,  p  =  »  and  (92)  becomes 


<A^>  -  <A>^   =  I  e^  <n^>  Cs^  A^ 
5  o 


(95) 


If  the  Initial  wavefront  is  a  point,  p  =  0  and  (92)  yields 


2       2   e    2    3  2 
<A  >  -  <A>  =  ~   <|i  >  Cs^  A 

15  o 


.[1] 


(9^) 


The  result  (95)  is  the  same  as  that  of  Krasilnlkov L J  vhile  (^k)    is  similar 
to,  but  not  identical  with,  that  of  P.G.  BergmannL  -^  . 


8.   The  reduced  vave  equation  in  a  random  medium 

— * 
Let  us  consider  the  solution  u(x)  of  the  following  problem 

Au  +  k  n  (x)u  =-6(x) 


11m   lx{ 


X   -►  00 


aixi 


-  iknu   =  0 


(95) 
(96) 
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Equation  (95)  is  the  reduced  wave  equation  with  a  source  term  corresponding 
to  a  source  of  unit  strength  located  at  the  origin  x  =  0.  Equation  (96)  is 
the  radiation  condition  which  asserts  that  the  solution  u  describes  outward- 
ly propagating  waves.   This  problem  is  well  posed  for  a  large  class  of  re- 
fractive Indices  n(x).   Let  us  assume  that  n(x)  has  the  form 

n(x)  =  1  +  €^I(1)  .  (97) 

Of  course  n,  (i  and  therefore  u  also  depend  upon  the  variable  w  which  ranges 
over  the  probability  space  Q,   but  we  shall  not  write  it  explicitly. 

It  seems  reasonable  to  suppose  that  u  can  be  represented  as  a  power  se- 
ries In  6 

u(x,€)  =  u^(x)  +  €u^(x)  +  e%(x)  +  0{e^)    .  (98) 

Upon  inserting   (97)    and   (98)    Into   (95)    ^<3^  {96),    and  equating  to  zero  the 
coefficient  of  each  power  of  e^   we  obtain 


2 
I     =  - 
o 

'Bu 


Au     +  k^u     =   -  6(x)  (99) 

o  o 


lim        |x|  ( -  iku    I     =0  (lOO) 


-    IKU      , 
O 


Au,    +  k^u,   =   -  2k^^u  (101) 

1  1  o 


(102) 


Au^  +  kV  =   -  k^n\   -   2k Vn  (1°5) 

2  2  o  1 

lim        1x1  ( -  iku^  -  iknuj     =  0   .  (10^) 
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The  solution  of  (99)  and  (lOO)  is 

iklxl 

%(^)  =  — — -  •  (105) 

^jt  1 X I 

To  solve  (lOl)  and  (102)  we  apply  Green's  theorem  to  the  interior  of  a  large 
sphere  and  assume  that  the  integral  over  the  sphere  tends  to  zero  as  its  ra- 
dius becomes  infinite.   Then  we  obtain 

2  p      lk|x-x'| 

\(^^=h     /  ^~^~^ n(x')ujx')d3^'  .  (106) 

^    |x-  x'  I 

In  a  similar  way  we  can  obtain  the  solution  u  (x)  of  (105)  and  (lO^),  which 

Is 

,2   p   ik|x-x'  I  p^  ~\ 

V^)=fe  /  ^  ^     ^  ^^(x')ujx')  +  2m(x')u^(x')  dS^'  .        (107) 


X-  X' 


Let  us  now  compute  the  mean  value  of  u  which  requires  the  calculation  of 
the  three  mean  values  <u  >,    <u^ >  and  <ap>.   From  (105)^  since  u  is  indepen- 
dent of  \i   and  therefore  of  w^  <u  >  =  u  .   From  (106),  since  u  is  linear  in  ^ 
and  since  we  assume  that  <|i>  =  0,  it  follows  that  <u^  >  =0.   To  compute  <Up> 
from  (107)  we  require  <n  >  and  <nu  >.  We  assume  that  <4i   (x)>  =  <^  >  is  a  con- 
stant^ and  we  compute  <^u^>by  using  (I06)  as  follows: 

,2   p   ik|x-x'| 

<n(x)u^(x)>=  1^   /  ^-^-^ <4i(x)M(x')>u^(x')d3^' 

^    |x  -  x' I 

,2^  2^  P   ik|x-x'| 

=  !L3L>     /  f N(lx-x'l)u  (x')dx' 

2«    J         Ix-x'l 

2  2    n   ikr 
=  ^-^    I     ^   N(r)u^(x+r)dr  .  (I08) 
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In  (108)  we  have  introduced  the  correlation  coefficient  N(|x-x'[)  = 
<^i(x)|i(x' )>/<n  >,   which  Is  assijmed  to  be  a  fimctlon  of  the  distance  |x-x'  | 
only. 

The  angular  integration  in  (108)  on  a  sphere  of  radius  r  can  be  per- 
formed explicitly  by  making  use  of  the  following  easily-proved  mean  value 
theorem  which  is  a  consequence  of  (99) 


1 


u  (x+r)ds 
o^ 


ujx) 


sin  kr 
kr 


0  <  r  <  |x| 


r  =  r 


u  (x)  sin  k|x| 


M 


r-  X 


kr 


r  >  X 


(109) 


When  (109)  is  inserted  into  (108)  it  becomes 


<^(x)u^(x)>=  u^(x)2kV 


e   N(r)sin  kr  dr 


i—  o 


-iklxl  .   ,  |-*|   r  21kr^,,  s^ 
+  e    I  'sm  k|'xl   /  e    N(r)±r 


=  -iu^(x)k<^  > 


/(e^^--l)N(r)-jr(e2i^^-l^|)-l)N(r)dr 


■-o 


(110) 

For  lx|  »  a^  which  holds  at  points  many  correlation  lengths  from  the  source, 
(110)  can  be  simplified  to 


<n(x)u^(x)>=  p<+i  >u^(x)  . 
Here  we  have  introduced  the  constant  3  which  is  defined  by 


(111) 


3  =  -ik   r  (e^^^''-l)N(r)dr  . 


(112) 
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Now  vhen  (ill)  is  used  in  it,  (10?)  yields 

iklx-x' ■ 


<U2(x)>=5-^^>{l+23)  /  ^—^-^ u^(x')d5J'  .  (115) 

J         Y  -  V  ' 


X  -  X 


Upon  collecting  our  results  we  have 

<u(x)>=  u  (x)  +  €  <u_(x)>+  0(e) 


2.  2  2         rv   iklx-x'  I 
,  u  (x)  +  ^^^"^  ^   (1+23)  / u  (x')d3'  +  0(£5).    (llif) 

Although  (11^)  gives  <u>  explicitly,  it  is  advantageous  to  rewrite  it  as  an 

2 
integral  equation  for  <u>.   Since  <u>  differs  from  u  by  0(€  ),  the  dlffer- 

2        2^ 
ence  between  e  <u>  and  e  u  is  0(e  )  and  we  may  write  (11^)  as 

2.  2  2         o   iklx-3c'  I 
<ai(x)>  =  u  (x)  +  ^\^  ^  (1+23)  /  ^-zr-Z —  ^^^'  )>  '^'  ^  °(^^)  •  (^^5) 

If  0(e    )    is  omitted,   this   equation  becomes   an  integral  equation  for  <u(x)>. 
It  is  equivalent  to  the   following  differential  equation  and  radiation  condi- 
tion for  <u> 

A<u>+  k^(n*)^<ii>=-6(x)  (116) 

lim        11^1  (^<^  -  ikn*u  )    =  0   .  (117) 

ri  VS|x|  J 

|x|    -»•   CO  ^11  ^ 

Here  the  effective   index  n*   is   defined  by 
(n*)^  =  1  +  €^  <^?>{\  +  23) 


00 


=  1  +  e^  <^?>   (  1  -  2ik   r(e^^^''-l)N(r)dr  J       .  (ll8) 


o 
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2 
Our  results  (ll6)  and  (11?)  shov  that  up  to  and  including  terms  of  order  e  , 

the  mean  value  of  u  propagates  as  If  in  a  medium  of  complex  refractive  index 

n*  given  by  (ll8).   This  is  the  main  result  of  this  section.   The  imaginary 

part  of  kn*  is  the  attenuation  coefficient  a  for  a  vave  in  this  medium.   From 

(ll8)  it  is  given  by 

00 

a  =  2e^  <n^>  k^  /  ^"^  "  '^°^   2kr)N(r)dr  .  (119) 

o 

This  result  agrees  with  that  obtained  previously  by  another  method.  L  J 

The  solution  of  (ll6)  and  (117)  is 

ikn*|x| 

<u(x>= .  (120) 

l..|xl 

By  comparing  (120)  with  (105)  we  see  that  <u(x)>  is  given  by  replacing  k  by 
kn*  in  u  (x) .   This  same  rule  applies  to  all  pairs  of  solutions  u  and  <u> 
of  problems  which  differ  from  (95),  (96)    or  (ll6),  (11?)  merely  in  the  source 
term.   From  (120)  we  have 

-alxl 

\<a(x)>\   = .  (121) 

W|x| 

The  method  which  we  have  used  to  derive  (ll6)  and  (11?)  can  also  be  used 

P  2  -+        ->■  2 

to  determine  <u  (x)>  and  from  it  the  variance  <u  (x)>  -  <u(x)>  .   The  quantity 

2 
<u>  is  often  called  the  coherent  wave  and  <u>  the  coherent  intensity  while 

the  variance  is  called  the  incoherent  intensity.   Since  it  should  be  clear  how 

2  2 

to  obtain  an  equation  for  <u  >  which  is  correct  throiigh  terms  in  €  ,  we  shall 

not  carry  out  this  calculation. 
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Caliromla  Division  Engineering  Library 

Dopartrrent  72-75,    Plant  A-1,    Bldg.   63-1 

Burbank,  California 

Attnt     N.  C.   Harnois 

The  Martin  Company 

P.  0.  Box  179 

Denver  1,  Coloraao 

Attn:   Mr.  Jack  McCormick 

The  Glenn  L.  Martin  Company 
Baltimore  3,  Maryland 
Attni   Engineering  Library 
Antenna  Design  Group 

Karvland  Electronic  Manufacturing  Corp. 
5009  Calvert  Road 
College  Park,  Maryland 
Attni  Mr.  H.  Warren  Cooper 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 


Plckard  and  Burns,  Inc. 

2ljO  Highland  Avenue 

Needhaa  %,  Mass. 

Attn:   Dr.  J.  T.  DeBettencourt 

Polytechnic  Research  and  Development 

Company,  Inc. 

202  Tillary  Street 

Brooklyn  1,  New  York 

Attn:  Technical  Library 

Radiation  Engineering  Laboratory 

Main  Street 
Maynard,  Mass, 
Attn:   Dr.  John  Ruze 

Radiation^  Inc. 

P.  0,  Drawer  37 

Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.L.  Cox 

Radio  Corp.  of  America 

RCA  Laboratories 

Rocky  Point,  New  York 

Attn:   P.  S,  Carter,  Lab.  Library 


Ryan  Aeronautical  Company 
Lindbergh  Field 
San  Diego  12,  Calirornl* 
Attni  Library  -  unclassified 

3age  Laboratories 
1?9  Linden  Street 
Wellesley  81,  Mass, 

Sanders  Associates 
95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  N.  R.  Wild,  Library 

Sandia  Corporation,  Sandla  Base 

P.O.  Box  5800,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  W.  Tumbull,  Engr.  Librarian 

Stanford  Research  Institute 

Menlo  Park,  California 

Attn:  Library,  Engineering  Division 


The  W,  L.  Maxscn  Corporation 

b60  West  "^bth  Street 

New  York,  N.  Y. 

Attn:     Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation 

Lambert  Saint-Louis  Municipal  Airport 

Box  516,   St.  Louis  3,   Missouri 

Attn:     R.  D.  Detrlch,   Engineering  Library 


Inc. 


McMillan  Laboratory 
Prownvllle  Avenue 
Ipswich,  Massachusetts 
Attn:   Security  Officer, 


Document  Room 


Melpar,    Inc.  (2) 

300Cp  Arlington  Boulevard 

Falls  Church,    Virginia 

Attn:     Engineering  Technical  Library 

Microwave  Develooment  Laboratory 
90  Broad  Street 
Babson  Park  57,   Massachusetts 
Attn:     N,  Tucker,   General  Manager 

Microwave  Ra  iiatlon  Company  Inc. 

19223  South  Hamilton  Street 

Gardena,   California 

Attn:     Mr.  Morris  J.  Ehrlich,    Prpsident 

Chance  Vought  Aircraft,    Inc. 
9311i  West  Jefferson  Street 
Dallas,   Texas 
Attni     Mr.   H.   S.  White,    Librarian 

Northrop  Aircraft,    Inc. 
Hawthorne,   California 

Attn:     Mr.  E.  A.  Freitas,   Library  Dept31U5 
1001  E.  Broadway 

Remlngt-on  Rand  Univ.   -  Division  of  Sperrv 

Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,    Pennsylvania 
Attn:     Mr.   John  F.   McCarthy 

R  and  D  Sales  and  Contracts 

North  American  Aviation,    Inc. 

1221ii  Lakewood  Boulevard 

Downey,   Califorr^la 

Attn:     Engineering  Library  b95-ll5 

North  American  Aviation,   Inc. 
los  Angeles   Intemati'^nal  Airport 
Los  Angeles  Ii5,    California 
Attn:     Engineering  Technical  Fiie 

Page  Comnunlcetions  Engineers,    Inc. 
710  Fourteenth  Street,   Northwest 
Washington  5,   D,  C. 
Attn:     Librarian 

Phllco  Corporation  Research  Division 

Branch  Library 

U700  Wlssachlckon  Avenue 

Philadelphia  biU,    Pa, 

Attnt     Mrs.  Dorothy  S.  Collins 


RCA  Laboratories 
David  Sarnoff  Research  Center 
Princeton,    New  Jersey 
Attn:     Miss  Fern  Cloak,    Librarian 
Research  Library 

Radio  Corporation  of  America 
Defense  Electronic  Prciucts 
Building  10,   Floor  7 
Camden  2,   New  Jersey 
Attn:     Mr,   Harold  J,   Schrader 

Staff  Engineer,  Organization 

of  Chief  Technical 

Adm'^nlstrator 

The  Ramo-Wooldridge  Corporation 
P.O,  Box  l45'453  Alrnort  Station 
Los  Angeles  Ii5,    California 
Attn:     Margaret  C,  Whitnah, 
Chief  Librarian 

Hoover  Microwave  Co. 
9592  Baltimore  Avenue 
College  Park,    Maryland 

Director,   USAF  Project  RAND 

Via:     Air  Force  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,    CallforrJ.a 

Rantec  Corporation 

Calabasas,   California 

Attn:     Grace  Keener,   Office  Manager 

Raytheon  Manufacturing  Company 
Miisslle  Systems  Division 
Bedford,    Mass. 
Attn:     Mr.   Irving  Goldstein 

Raytheon  Manufacturing  Company 
Way! and  Laboratory,    State  Road 
Wayland,   Mass. 
Attn:     Mr.  Robert  Borts 

Raytheon  Manufacturing  Comr«ny 
Wayland  Laboratory 
Wayland,   Mass. 

Attn:     Miss  Alice  G,  Anderson, 
Librarian 

Republic  Aviation  Corporation 
Farmingdale,    Long   Island,   N.  Y. 
Attn:     Engineering  Library 

Thru:     Air  Force  ^lant  Representative 
Republic  Aviation  Corp. 
Famlngdale,   ^^ng  Island,  N»Y, 

Rheem  Manufacturing  Conpany 
9236  East  Hall  Road 
Downey,    California 
Attn:     J.  C.   Joerger 

Trana-Tech,    Inc. 
P.  0.  ^ox  31*6 
Frederick,   Maryland 


Sylvania  Electric  Products,    Inc. 
100  First  Avenue 
Walthara  51i,   Mass, 

Attn:     Charles  A,  Thornhill,   Report  Librarian 
Waltham  Laboratories  Library 

Systems  Laboratories  Corporation 
llifl52   Ventura  Boulevard 
Sherman  Oaks,   California 
Attn:     Donald  L.  Margerum 

TRG,  Inc. 

17  Union  Square  West 

New  York  3,  N.  Y, 

Attn:     M,   L,   Henderson,    Librarian 

A.   S,   ThOTias,    Inc. 

I6l  Devonshire  Street 

Boston  10,   Mass. 

Attn:  A.  S.  Thomas,  President 

Bell  Telephone  Laboratories 
Miirray  Hill 
New  Jersey 

Chu  Associates 

P.  0.  Box  3P7 
Whitcomb  Avenue 
Littleton,  Mass. 

Microwave  Associates,  Inc. 
Burlington,  Mass. 

Raytheon  Manufacturing  Company 
Missil**  Division 
Hartwell  Road 
Bedford,  Mass. 

Radio  Corporation  of  America 
Aviation  Systems  L-aboratory 
225  Crescent  Street 
Waltham,  Mass, 

Lockheed  Aircraft  Corporation 
Missile  Systems  Division  Research  Library 
Box  50U,  Sunnyvale,  California 
Attn:  Miss  Eva  Lou  Robertson, 
Chief  Librarian 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:   Dr.  W.  C.  Hoffman 

Ccvqmander 

AF  Office  of  Scientific  Research 
Air  Research  and  Development  Comand 
■"it^  Street  and  Constitution  Avenue 
Washington,  D.  C. 
Attn:  Mr.  Ottlng,  SRY 

Westinghouse  Electric  Corp. 
Electronics  Division 
Friendship  Int'l  Airport  Box  7li6 
Baltimore  3,  Maryland 
Attn;   Engineering  Library 
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Wheeler  laboratories.    Inc. 
122  Cotter  Mill  Road 
Great  Neck,   New  York 
Attn:     Mr.  Harold  A.  Wheeler 

Zenith  Plastics  Co. 
Box  91 

Gardena,   California 
Attni     Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute 
of  the  Dnlversity  of  Alaska 
College 
Alaska 

University  of  California 

Berkeley  U,   California 

Attnt     Dr.  Saaiuel  Silver, 

Prof.  Engineering  Science 
Division  of  Elec.  Eng.   Electronics 
Research  Lab. 

UniTersity  of  California 
Electronics  Research  Lab. 
332  Cory  Hall 
Berkeley  h,  Califomia 
Attn:     J.  R.  Whinnery 

Califomia  Institute  of  Technology 
Jet  Propulsion  Laboratory 
UROO  Oak  Orove  Drive 
Pasadena,  California 
Attn:     Mr.    I.   E.   Newlan 

California  Institute  of  Technology 
1201  E.  Califomia  Street 
Pasadena,   California 
Attn:     Dr.  C.   Papas 

Carnegie  Institute  of  Technology, 
Schenley  Park 

Pittsburgh  13,    Pennsylvania 
Attn:     Prof.  A.   E.   Heins 

Crnell  University 

School  of  Electrical  Engineering 

Ithaca,   New  York 

Attnt     Prof.  G.   C.   Dalnan 

I'niversity  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,   Florida 

Attn:      Prof.   M.   H.   Latour,    Library 

Lih'-ary 

Georgia  Institute  of  Technology 

Engineering  Experiment  Station 

Atlanta,   Georgia 

Attn:     Mrs.   J.H.   Crosland,    Librarian 

Harvard  University 

Technical  Reports  Collection 

Gordon  McKay  Library,   303A  Pierce  Hall 

Oxford  Street,   Cambridge  38,  Mass. 

Attn:     Mrs.  E.L.  Hufschmidt,   Librarian 

Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39,   Mass. 
Attn:     Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,    Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,    Illinois 

Attn:     Dr.   P.  E.   Moyes,   Department  of 
Electrical   Engineering 

Tlia- Johns  Hopkins  University 
Homewood  Campus 
Department  of   Physics 
Baltimore  18,  Maryland 
Attn:     Dr.  Donald  E.  Kerr 

Sandia  Corporation 
Attn:     Organization  1[|23 
Sandia  Base 
Albuouerqoe,   New  Mexico 


Applied  Physics  Laboratory 
The  JohflSHopkins  University 
8621  Georgia  Avenue 
Silver  Spring,   Maryland 
Attn:     Mr.   George  L.  Seielstad 

Massachusetts  Institute  of  Technology 

Research  Laboratory  of  Electronics 

Room  20B-221 

Cambridge  39,  Massachusetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology 

Lincoln  Laboratory 

P.  0.  Box  73 

Lexington  73,  Mass. 

Attn:  Document  Room  A-229 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J,  A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 
Attn:  Prof.  K.  M.  Siegel 
912  N.  Main  St., 
Ann  Arbor,  Michigan 

Unlversitv  of  Michigan 
Engineering  Research  Institute 
Willow  Run  Laboratories 
Willow  Run  Airport 
Ypsilantl,  Michigan 
Attn:  Librarian 

University  of  Minnesota 

Minneapolis  Id,  Minnesota 

Attn:  Mr.  Robert  H.  Stumm,  Library 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:   R.  E.  Bean 

Ohio  State  University  Research  Found. 
Ohio  State  University 
Columbus  10,  Ohio 
Attn:  Dr.  T.E.  Tice 

Dept.  of  Elec.  Engineering 

The  UniTersity  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attni   Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn,  New  York 

Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
5?  Johnson  Street 
Brooklyn.  New  York 
Attn:  Mr.  A.  E.  Laemmel 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  New  York 
Attn:  Dr.  C.  S.  Grove,  Jr. 

Director  of  Engineering  Research 

The  University  of  Texas 
Elec.  Engineering  Research  Laboratory 
P.  0.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R.  Gerhardt 
Assistant  Director 

The  University  of  Texas 

Defense  Research  Laboratory 

Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 

Department  of  Electrical  Engineering 

Toronto,  Canada 

Attn:   Prof,  G.  Sinclair 


Lowell  Technological  Institute 
Research  Foundation 
P.  0.  Box  709,  Lowell,  Mass. 
Attn:  Dr.  Charles  R.  Minglns 

University  of  Washington 

Department  of  Electrical  Engineering 

Seattle  5,  Washington 

Attn:  G,  Held,  Associate  Professor 

Stanford  University 
Stanford,  Califomia 
Attn:  Dr.  Chodorow 

Microwave  Laboratory 

Physical  Science  Laboratory 

New  Mexico  College  of  Agricultiire 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr.  H.  W.  Haas 

Brown  University 
Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr.  C.  M.  Angulo 

Case  Institute  of  Technology 

Cleveland,  Ohio 

Attn:  Prof.  S.  Seeley 

Columbia  University 

Department  of  Electrical  Engineering 

Momlngside  Heights 

New  York,  N.  Y. 

Attn:  Dr.  Schlesinger 

McGill  University 

Montreal,  Canada 

Attn:  Prof.  0.  A.  Woonton 

Director,  The  Eaton  Electronics 

Research  Lab. 

Purdue  University 

Department  of  Electrical  Engineering 

Lafayette,  Indiana 

Attn:  Dr.  Schulta 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania 

Institute  of  Cooperative  Research 

3bOO  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Dept,  of  Electrical  Engineering 

University  of  Tennessee 
Ferris  Hall 

W.  Cumberland  Avenue 

Knoxville  16,  Tennessee 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

University  of  Seattle 

Department  of  Electrical  Engineering 

Seattle,  Washington 

Attn:  Dr.  D.  K.  Reynolds 

Wayne  University 

Detroit,  Michigan 

Attn:  Prof.  A.  F.  Stevenson 

Electronics  Research  Laboratory 
Illinois  Institute  of  Technology 
3300  So.  Federal  Street 
Chicago  16,  Illinois 
Attn:  Dr.  Lester  C,  Peach 
Research  Engineer 

Advisory  Group  on  Electronic  Parts 

Room  103 

Moore  School  Building 

200  South  33rd  Street 

Philadelphia  U,  Pennsylvania 
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Ionosphere  Research  Laboratory 

Pennsvlvanla  State  College 

State  College,    Pennsylvania 

ATTMi     Professor  A.   H.  Waynlck,   Director 

Institute  of  Mathematical   Sciences 

25  Wayerly  Place 

New  Tork  3,   New  York 

ATTN  I     Librarian  ^3' 

Electronics  Division 
Rand  Corporation 
1700  Main  Street 
Santa  Monica,   California 
ATTNt     Dr.   Robert  Kalaba 

National  Bureau  Df  Standards 

Washington,   D.  C. 

ATTNt     Dr.  W.   K.  Saunders 

Applied  Mathenatlcs  and  Statist  irs  La*^, 

Stanford  University 

Stanford,     ttlifomia 

ATTN;     Dr.  Albert  H.  Bowker 

Department  of  Physics  and  A3tr:nomy 
Michigan  State  College 
East  Lansing,   Michigan 
ATTN:     Dr.  A.  Leltner 

University  of  Tennessee 
Knoxvllle,  Tennessee 
ATTN!     Dr.  Fred  A.  Flcken 

Lebanon  Valley  College 
Annvllle,  Pennsylvania 
ATTN:     Professor  B.H.  Bissinger 

General  Atomic 

P.  0.  Box  608 

San  Dleeo  12,  California 

Attn:      fV.   Edward  Gerjuoy 

Department  of  Physics 
Amherst  College 
Anherst,   Mass. 
ATTN:     Dr.  Arnold  Arons 

California  Institute  of     Technology 
1201  E.  California  Street 
Pasadena,   California 
ATTN!     Dr.   A.  Erdelyl 

Mathematics  Department 

Stanford  Hniversity 
Stanford,  California 
ATTN:     Dr.   Harold  Levins 

University  of  Minnesota 

Mlnnpapolls  I'l,   Minnesota 

ATTNt     Professor  Paul  C.   Rosenbloom 

Department  of  Mathematics 

Stanford  University 

Stanford,   California 

ATTN:     Professor  Bernard  Epstein 

Applied  Physics  Laboratory 
The  Joh* Hopkins  University 
8621  Georgia  Avenue 
Silver  Spring.   Maryland 
ATTN:     Dr.  B.  S.  Gourary 

(2)Exchange  and  Gift  Division 
The  Library  of  Congress 
Washington  25,    D.  C. 

Electrical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,    Mass. 
ATTNt     Dr.   L.   J.  Chu 

Nuclear  Developnent  Associates,    Inc. 

5  New  Street 

White  Plains,   New  York 

ATTN:     Library 

California  Institute  of  Technology 
Electrical  Engineering 
Pasadena,   California 
ATTNt     Dr.   Zohrab  A.  Kaprlellan 


Dr.   Rodman  Doll 
311  W.  Cross  Street 
Tpsllantl,  Hlchigan 

Callfoimla  Inst,   of  Technology 
Pasadena,    Callfon^ia 
ATTN:     Mr.   Calvin  Wilcox 

Mr.  Robert  Brockhurst 

Woods  Hole  Oceanographlc  Institute 

Woods  Hole,   Mass. 

National  Bureau  of  Standards 
Boulder,   Colorado 
ATTN:     Dr.   R.  Oallet 

Dr.   Solomon  L.  Schwebel 

36fi9  Louis  Road 

Palo  Alto,     Callfomia 

University  of  Minnesota 
The  University  of  Library 
Minneapolis  lU,   Minnesota 
ATTNt     Exchange  Division 

Depai^ment  of  Mathematics 
University  of  Callfomia 
Berkelev,   California 
ATTM:      Professor  Bernard  Friedman 

Lincoln  Laboratory 

Massachusetts  Institute  of  Technology 

P.  0.  Box  73 

Lexington  73,   Massachusetts 

ATTN:     Dr.  Shou  Chin  Wang,    Room  C-351 

Melpar,  Inc., 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

ATTN:  Mr.  K.  S.  Kelleher,  Section  Head 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  u.  Hanscom  Field 

Bedford,  Mass. 

ATTN:  Hr.  Francis  J.  Zuckar,  CRHD 

Hq.  Air  Force  Cambridge  Research  Center 

Laurence  G.  Hanscom  Field 

Bedford,  Mass. 

ATTN:   Dr.  Philip  Newman,  CRRK 

Mr.  N.  C.  Cerson 

Trapelo  Road 

South  Lincoln,  Mass. 

Dr.  Richard  B.  Ban-ar 

Systems  Development  Corp. 
2(:00  Colorado  Avenue 
Santa  Monica,  California 

Columbia  University  Hudson  Laboratories 
P.O.  Box  239 

1U5  Palisade  Street,  Dobbs  Ferry,  N.  Y. 
ATTNt  Dr.  N.  W.  Johnson 

Institute  of  Fluid  Dynamics 
and  Applied  Mathematics 
University  of  Maryland 
College  Park,  Maryland 
ATTNt  Dr.  Elliott  Montroll 

Department  of  Electrical  Engineering 

Washington  University 

Saint  Louis  5,  Mo. 

ATTN:     Professor  J.   Van  Bladel 

DepaiH-raent  of  the  Navy 

Office  of  Naval  Re'iearch  Branch  Office 

1030  E.  Green  Street 

Pasadena  1,   Callfomia 

Brandels  ^Tniverslty 
Waltham,   Mass. 
ATTN:     Library 


General  Electric  Company 
Mictowave  Laboratory 
Electi^nics  Division 
Stanford  Industrial  Park 
Palo  Alto,   California 
ATTN:     Library 


Snyth  Research  Associates 
3555  Aei^)  Court 
San  Diego  3,  Callfomia 
ATTNt     Dr.   John  B.  Smyth 

Electrical  Engineering 
California   Institute  of  Technology 
Pasadena,   California 
«TTS:     Drfleorges  G.  Weill 

Naval  Research  Laboratory 

Washington  25,    D.  C. 

ATTN:     Henry  J.  Passerlnl,   Code  5278A 

Dr.  George  Kear 
5  Culver  Court 
Oi^nda,   California 

Brooklyn  Polytechnic 

85  Livingston  Street 

Brooklyn,    New  York 

ATTNt     Dr.  Nathan  Marcuvlti 

Department  of  Electrical  Engineering 
Brooklyn  Polytechnic 
P?  Livingston  Street 
Brooklyn,   New  York 
ATTN:     Dr.   Jerry  Shmoys 

Department  of  Mathematics 
University  of  New  Mexico 
Albuquerque,    New  Mexico 
ATTN:     Dr.    I.  tolodner 

Mathematics  Department 
Polytechnic   Institute  of  Brooklyn 
Johnson  and  Jay  Street 
Brooklyn,   lew  York 
ATTN:     Dr.  Harry  Hochstadt 

Ballistics  Research  Laboratory 
Aberdeen  Pi^vlng  Grounds 
Aberdeen,   Maryland 
ATTN:     Dr.   Pullen  Keats 

Dr.  Lester  Kraus 
):935  Whitehaven  Way 
San  Diego,   California 

University  of  Minnesota 
Institute  of  Technology 
Minneapolis,   Minnesota 
Attn:     Dean  Athelston  Spilhaus 

Ohio  State  University 
Columbus,   Ohio 
Attn:     Prof.  C.  T.  Tai 
Depai^ment  of  Electrical  Eng. 

Naval  Research  Laboratories 

Washington  25,   D.  C. 

Attn:     W.   S.   Ament,   Coda  5271 

Naval  Research  Laboratory 
Washington  25.    D.  C, 
Attn:     Dr.  Leslie  G.  McCracken,    Jr. 
Code  3933A 

Office  of  Naval  Research 

Department  of  the  Navy 

Attn:     Geophysics  Branch,   Code  Iil6 

Washington  25,   D.  C. 

Office  of  Chief  Signal  Officer 
Signal  Plans  and  Operations  Division 
Attn:     SIOOL-2.    Room  20 
Con.   Liaison  Br.,    Radio  Prop.   Sect. 
The  Pentagon,   Washington  25,    D.   C. 

Defence  Research  Member 
Canadian  Joint  Staff 
2001  Connecticut  Street 
Washin^on,    D.   C. 

Central  Radio  Prop.   Lab. 
National  Bureau  of  Standards 
Attn:     Technical  Reports  Library 
Boulder,   Colorado 

D.  S.  Weather  Bureau 
U.  S.  Depairtnent  of  Comnerce 
Washlnpton  25,   D.   C. 
Attn:     Dr.   Harry  Wexler 
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Federal  Communications  Commission 

Washington  25,    D.   C. 

Attn:     Mrs.    Barbara  C.   Grimes,    Librarian 

Upper  Atmosphere  Research  Section 
Central  Radio  Propagation  Laboratory 
National  Bureau  of  Stareiards 
Boulder,   Colorado 

Argonne  National  Laboratory 

P.O.  Box  299 

Lemont ,  Illinois 

Attn:  Dr.  Hoylande  D.  Young 

Bell  Telephone  Labs. 

Murray  Hill,  New  Jersey 

Attn:  Dr.  S.  0.  Rice,  33  -  203 

Carnegie  Institute  of  Washington 
Dept.  of  Terrestrial  Magnetism 
52I1I  Broad  Branch  Road,  N.  W. 
Washington  1?,  D.  C. 
Attn:  Library 

Georgia  Tech  Research  Institute 
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